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ABSTRACT 

We decompose the rotation curves of 34 bright spiral galaxies into baryonic and dark 
matter components. Stellar mass profiles are created by applying color-M/L relations 
to near-infrared and optical photometry. We find that the radial profile of the baryonic- 
to-dark- matter ratio is self-similar for all galaxies, when scaled to the radius where the 
■ contribution of the baryonic mass to the rotation curve equals that of the dark matter 

(R x ). We argue that this is due to the quasi-exponential nature of disks and rotation 
curves that are nearly flat after an initial rise. The radius Rx is found to correlate 

C\| ' most strongly with baryonic rotation speed, such that galaxies with Rx measurements 

that lie further out in their disks rotate faster. This quantity also correlates very 
strongly with stellar mass, Hubble type, and observed rotation speed; S-band central 
surface brightness is less related to Rx than these other galaxy properties. Most of the 
O galaxies in our sample appear to be close to maximal disk. For these galaxies, we find 

that maximum observed rotation speeds are tightly correlated with maximum rotation 
speeds predicted from the baryon distributions, such that one can create a Tully-Fisher 
relation based on surface photometry and redshifts alone. Finally, we compare our 
data to the NFW parameterization for dark matter profiles with and without including 
adiabatic contraction as it is most commonly implemented. Fits are generally poor, and 

^ ■ all but 2 galaxies are better fit if adiabatic contraction is not performed. In order to 

have better fits, and especially to accommodate adiabatic contraction, baryons would 
need to contribute very little to the total mass in the inner parts of galaxies, seemingly 
in contrast with other observational constraints. 

Subject headings: galaxies: fundamental parameters - galaxies: general - galaxies: halos - 
galaxies: kinematics and dynamics - galaxies: spiral - galaxies: stellar content 
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heijen 1997; de Block & McGaugh 1998). In- 
terestingly, the main uncertainty in the radial 
density distribution of dark matter, as derived 
from rotation curves, stems from the poorly 
known stellar mass distribution (e.g., Verhei- 
jen 1997). 

In this paper, we do not rely on the typi- 
cal maximal disk assumption (van Albada et 
al. 1985; van Albada Sz Sancisi 1986) to de- 
rive stellar mass profiles, but use color- mass- 
to-light ratio (M/L) relations first given by 
Bell & de Jong (2001) and later updated by 
Bell et al. (2003). These color- Af/L relations 
were determined from the analysis of spectro- 
photometric spiral galaxy evolution models 
and give an upper limit to the baryonic mass 
present in spirals. They are most power- 
ful when applied to data in the near-infrared 
since this wavelength regime traces the older 
stellar populations that contain most of the 
mass, while avoiding much of the effects of 
dust (extinction at the near-infrared i^T-band 
is only about 10% of that at the optical B- 
band; Martin & Whittet 1990). The color- 
M/L relations allow us to investigate radial 
variation in stellar M/L, and provide us with 
a consistent way of scaling stellar M/L from 
one galaxy to the next. With the distributions 
of stellar mass in galaxies reasonably well de- 
termined, we can begin to investigate those of 
the dark matter. 

The maximum disk assumption (Sackett 
1997) can be examined with the color-M/L 
relations. We can obtain an upper limit to 
the number of galaxies that have maximal 
disks, and also compare the shapes of rota- 
tion curves derived from baryonic mass dis- 
tributions in galaxies to the observed rota- 
tion curves as in Palunas & Williams (2000). 
There are many pieces of evidence in the lit- 
erature for maximal disks; we mention a few 
noteworthy ones here. Possibly the strongest 
evidence comes from the work of Weiner, Sell- 
wood, & Williams (2001) who used fluid- 



dynamical models of gas flow to model the 
two-dimensional velocity field of a barred 
galaxy. They find that such models require 
80 - 100% of the maximal disk M/L[ value. 
Palunas & Williams (2000) modeled the bary- 
onic mass distributions of 74 spirals and found 
that a mass-follows-light model could repro- 
duce the overall structure of the optical rota- 
tion curves in the majority of galaxies. These 
authors found that 75% of the galaxies in 
their sample have a rotation curve out to -R23.5 
that is entirely accounted for by baryons. The 
mass models for 20% of their sample fail be- 
cause of nonaxysymmetric structures. For a 
more inclusive discussion on maximal disks, 
see Palunas & Williams (2000). 

It has been known for nearly 20 years that 
the radius where dark matter begins to con- 
tribute to the rotation curve of a galaxy is 
smaller for low surface brightness galaxies 
than it is for brighter galaxies (e.g., Persic 
& Salucci 1988, 1990; Broeils 1992; Persic, 
Salucci, & Stel 1996; de Blok & McGaugh 
1997). A number of authors have demon- 
strated this by evaluating either the mass dis- 
crepancy of galaxy disks (ratio of total to 
baryonic mass) by adopting a stellar M/L, or 
a similar quantity at a chosen radius. For ex- 
ample, Persic, Salucci, & Stel (1996) find that 
spiral galaxies with V = 100 km s _1 have 
> 75% of their mass in dark matter within 
Ropt = 3.2/i where h is the disk scale-length, 
and that galaxies with V = 150 km/s have 
only > 40%. McGaugh & de Blok (1998) 
find that the total M/Lb evaluated at Ah is 
smaller for galaxies of higher surface bright- 
ness and brighter magnitude, and does not 
correlate with h. These authors also find that 
the radius where dark matter begins to dom- 
inate over the luminous mass is greater for 
higher surface brightness galaxies. Zavala et 
al. (2003) find that the ratio of the velocity 
due to the baryons to the total maximum ve- 
locity depends mainly on disk central surface 



2 



density such that denser galaxies have larger 
ratios. These authors also find that the mass 
discrepancy evaluated at 2.2h and 5h depends 
on central surface brightness, but puzzlingly 
does not depend on h, baryonic mass, or B- 
band luminosity. In addition, Zavala et al. 
(2003) find that M /Lb anti-correlates signifi- 
cantly with S-band central surface brightness 
and does not correlate with 5-band luminos- 
ity, in disagreement with McGaugh k de Blok 
(1998) and the general qualitative trend in 
the literature. Most recently, Pizagno et al. 
(2005) evaluate the total M/Lj at 2.2/ij and 
find median values of 2.4 and 4.4 for galaxies 
with stellar masses greater than 1O 1O M and 
4.4 between 1O 9 -1O 1O M , respectively. These 
authors are all generally in qualitative agree- 
ment, but there are discrepancies, and we will 
discuss their origin in the body of this paper. 

A combination of resolved Ha or CO and 
H I velocity data is necessary to determine the 
dark matter distribution of a galaxy from its 
center to beyond ~ 5h. Ha data are essential 
to measure the steep rise of rotation curves 
in the inner ~ 2h of galaxies (e.g., Palunas 
k Williams 2000), whereas H I data are key 
to probe beyond where Ha can be measured 
(e.g., de Blok k McGaugh 1997). Since ro- 
tation curves in the optical portion of galax- 
ies are usually affected by non-axisymmetric 
structures that do not trace the main galac- 
tic potential, such as bars and spiral arms, 
the distributions of dark matter derived from 
Ha data can be quite variable. Two examples 
of studies probing dark matter distributions 
that have used both types of rotation curves 
in a direct manner, independent of prior pa- 
rameterizations for the dark matter, are Per- 
sic, Salucci, k Stel (1996) and McGaugh k 
de Blok (1998). Persic, Salucci, k Stel (1996) 
coadd 1100 rotation curves and find that they 
could be determined by a single galactic pa- 
rameter (e.g., luminosity). This formulation, 
however, is refuted by a number of authors 



(e.g., Verheijen 1997; Bosma 1981). McGaugh 
k de Blok (1998) find a regularity between 
the mass discrepancy and acceleration of a 
galaxy as a function of radius. This follows 
directly from the Tully-Fisher relation, prob- 
ably the best observational example of self- 
similarity in the dark matter component of 
galaxies. The mass discrepancy-acceleration 
relation shows that knowledge of the baryonic 
mass distribution of a galaxy allows one to 
calculate its dark matter distribution. Such a 
relation steps beyond the usual Tully-Fisher 
relation, which predicts the dark matter con- 
tent of a galaxy at a particular radius given 
its baryonic mass at that radius, to one that 
predicts radial distributions. 

Theories and simulations of galaxy forma- 
tion should be able to reproduce the distribu- 
tions of dark matter observed in galaxies. Un- 
til now, comparisons to predictions for dark 
matter distributions have consisted of fits of 
functional formulations for dark matter halos 
to rotation curve data, with stellar M/L as 
a free parameter. Now that we have a han- 
dle on the stellar M/L in galaxies with the 
color-M/L relations, we re-examine these fits. 
In particular, we test the main incarnation of 
the density distribution of dark matter halos 
in N-body simulations of a ACDM universe: 
the formulation of Navarro, Frenk, k White 
(1996) (NFW), which is the simplest and most 
popular analytical description. This formula- 
tion predicts that dark matter is significant 
down to the very inner radii of galaxies, in 
contradiction to maximal disks. Moreover, 
dark matter halos are expected to contract 
due to baryon collapse. As this contraction 
is normally implemented (Blumenthal et al. 
1986), even more dark matter is caused to 
move toward the centers of galaxies, making 
the situation worse (e.g., McGaugh k de Blok 
1998). We perform fits to the NFW formula- 
tion with and without including contraction 
as it is normally implemented. The main as- 
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sumption behind these fits, namely that the 
NFW halos can be fit to data that subtends 
only their very inner parts (< 0.1i?„j r j o /), is 
problematic. Furthermore, the simulations 
that produce the NFW functions have diffi- 
culty forming disk galaxies once hydrodynam- 
ics is included (e.g., Abadi et al. 2003), sig- 
naling possible unknown interactions between 
baryons and dark matter. To address these 
concerns, we turn the problem of fitting dark 
matter halos to data around and present a 
simple universal form for dark matter profiles 
in terms of baryonic mass profiles. 

This paper is organized as follows: In §2 
we briefly discuss our galaxy sample. Radial 
baryonic mass distributions and their rotation 
curves are derived in §3. In §4 these baryonic 
rotation curves are compared with rotation 
curves from the literature to derive rotation 
curves due to dark matter. We discuss max- 
imal disks in §5. In §6 we define quantities 
to describe the baryonic and dark matter dis- 
tributions and show how they correlate with 
general galaxy properties. We investigate the 
radial behavior of dark matter in §7, and mass 
profiles are fit with the NFW profile in §8. In 
§9, we summarize our conclusions. Through- 
out this paper we adopt a Hubble constant 
(H Q ) of 70 km s _1 Mpc -1 . When distance- 
dependent quantities have been derived from 
the literature, we have reduced them to this 
value of H Q and always quote the converted 
values. 

2. The Galaxy Sample 

Our data consist of surface brightness pro- 
files, physical parameters, and rotation curves 
for 34 bright spiral galaxies. These galaxies 
have inclinations in the range ~ 30-65 degrees 
in order to reduce the effects of dust, while 
still being able to obtain accurate kinematical 
information. Surface brightness profiles for 30 
galaxies were presented in Kassin, de Jong, 
& Pogge (2005) (hereafter, Paper I); those for 



the remaining 4 galaxies can be found in Ver- 
heijen (1997). More details on the sample 
selection are given in Paper I. The profiles 
were calculated in elliptical annuli of increas- 
ing distance from the centers of the galaxies 
and have been corrected for Galactic extinc- 
tion. The galaxies NGC 1090, NGC 2841, and 
NGC 3198, which have images from the Sloan 
Digital Sky Survey (SDSS) Second Data Re- 
lease (Abazajian et al. 2004), have images in 
which almost half of the galaxy was off of 
the detector. This is also the case, but to 
a much lesser extent for NGC 3521. This ef- 
fect can be observed in the g-band images in 
Figure 1 of Paper I; these galaxies are flagged 
in the following analysis. All the analyses in 
this paper have been performed for both the 
approaching and receding sides of a handful 
of galaxies with complete imaging. No sig- 
nificant change (to within the zero-point un- 
certainties) was found between the results for 
each side of these galaxies. Physical param- 
eters (i.e., Hubble type, distance, R25, inte- 
grated magnitudes) and bulge-disk decompo- 
sitions can also be found in Paper I. 

3. Baryonic Matter 

In this section, we derive the baryonic mass 
surface density profiles of the galaxies, and 
compute the component of the galaxies' ob- 
served rotation that is due to the baryons. 
These components are called "baryonic rota- 
tion curves," to distinguish them from ob- 
served rotation curves. In §4, we will intro- 
duce the term "dark matter rotation curves," 
to denote the contribution from dark matter 
to the observed rotation curves. 

3.1. Radial Baryonic Surface Mass- 
Density Distributions 

To determine a galaxy's radial stellar sur- 
face mass-density distribution, we apply a 
color-M/L relation to its surface brightness 
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Fig. 1. — Color-M/L relation for B — R color and 
(M/L), at K from Bell ct al. (2003) (solid line). 
The galaxies in the sample are plotted as filled 
circles, and a reddening vector is plotted for for 
correction to face-on for a Milky Way-type galaxy 
viewed at an inclination of 80° with the Tully et 
al. 1998 formalism. 

profiles. These relations were derived from 
spectrophotometric spiral galaxy evolution 
models in Bell Sz de Jong (2001), and were 
updated in Bell et al. (2003). The color-M/L 
relations show a relatively tight correlation 
(~ 0.1 dex spread for the color and metallic- 
ity range where our galaxies lie) between the 
optical color of a galaxy and its stellar mass- 
to-light ratio, (M/L)*. Color-M/L relations 
are most useful when applied to (M/L)* in the 
near-infrared since they are the least affected 
by dust obscuration, as discussed in Bell & 
de Jong (2001). Specifically, we choose to use 
the relation between B — R color and (M/L)* 
in the if-band. This relation is composed of 
the optical color with the largest wavelength 
baseline and the reddest near- infrared band. 
It is reproduced in Figure 1. If imaging is un- 
available at R, we use the relation for B — V, 
and if imaging is unavailable at K we use H. 

For those galaxies without a significant 



bulge contribution (such that inclusion of a 
bulge component does not change the rota- 
tion curve due to baryons beyond the uncer- 
tainties of the color-M/L relations), we ap- 
ply a color-M/L relation directly to the az- 
imuthally averaged radial B — R color profiles 
to derive radial (M/L)* profiles at K. The 
Lk profile for each galaxy is then multiplied 
by the galaxy's (M/L)* profile to derive a ra- 
dial stellar surface mass-density profile. For 
these galaxies, using a radial B — R color pro- 
file is consistent with using an aperture B — R 
color since this color-M/L relation has a shal- 
low slope. For those galaxies with a signifi- 
cant bulge component, a bulge-disk decompo- 
sition is performed. A characteristic B — R 
color is adopted for each component based on 
the average colors of the bulge and disk, and 
the Lk profiles are multiplied by the resulting 
(M/L)*. Extinction corrections are discussed 
in §3.4. We also extend the exponential disk 
surface brightness profiles until approximately 
oo (defined here as r = 10000" ). 1 

When available from the literature, we use 
radial H I 21 cm measurements to determine 
the contribution of interstellar gas to the ra- 
dial baryonic surface mass-densities. There 
are gas measurements available for NGC 1090 
from Gentile et al. (2004), NGC 3198 from 
Begeman (1987), and NGC 3949, NGC 3953, 
and NGC 3992 from Verheijen (1997). The 
neutral gas component is included by scal- 
ing the H I surface mass-density by a fac- 
tor of 1.32 to account for the abundance 
of helium. For those galaxies without gas 
mass measurements, we assume that gas does 
not contribute to its baryonic surface mass- 
density. This is probably not a bad assump- 
tion for high surface brightness galaxies since 
a galaxy with a Hubble type between Sa and 
Scd has a gas mass fraction that is typically 

1 Extended galactic stellar disks have been discovered 
for some galaxies via very deep imaging (see Ibata et 
al. 2005, and references within). 
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Fig. 2 — Rotation curve for NGC 1090 due to the 
stellar mass component (solid line) is compared 
with the rotation curve due to stars and H I gas 
from Gentile ct al. 2004 (dashed line). The ob- 
served Ha and H I rotation curves are plotted as 
open triangles and circles, respectively. The up- 
per and lower bounds for the stellar mass rotation 
curve (dotted lines) are due to the 0.1 dcx uncer- 
tainty in the color-M/L relations. 

M HI /M gas+Stars ps 0.03 (Roberts & Haynes 
1994), and hence does not greatly affect its 
baryonic mass. The average gas mass frac- 
tion for the 5 galaxies in our sample that 
have gas mass measurements is 0.04, consis- 
tent with this measurement. Figure 2 demon- 
strates the effect of including the interstel- 
lar gas component in the baryonic rotation 
curve of NGC 1090, which is generally within 
the uncertainty of the stellar mass rotation 
curve. Therefore, by not including interstel- 
lar medium contributions in the analysis, only 
a minor systematic effect is introduced. 

3.2. Baryonic Rotation Curves 

A baryonic rotation curve is calculated for 
each galaxy from its baryonic surface mass- 
density profile and is plotted in Figure 3. The 
rotation curve calculation is done with the 



rotmod task in the GIPSY software package 
(van der Hulst et al. 1992), which calculates 
rotation curves for galaxies composed of a 
truncated exponential disk (Casertano 1983) 
and a spherical bulge when applicable. Ex- 
ponential disks are assumed to have a scale- 
height of 0.3 kpc, which is typical of bright 
spirals. The difference between using a differ- 
ent scale-height for each galaxy of O.lhm kpc, 
as in Sparke & Gallagher (2000), and a global 
value of 0.3 kpc is negligible. The results also 
remain unchanged if we used the relation be- 
tween the central surface brightness and the 
ratio of vertical scale-height to vertical scale- 
length given by Bizyaev & Mitronova (2002). 

3.3. Uncertainties in Baryonic Rota- 
tion Curves 

The largest source of uncertainty for the 
baryonic rotation curves lies in the normal- 
ization of the color-M/L relation, which is 
mainly determined by the stellar IMF at low- 
mass. Since the faint end of the IMF is rel- 
atively unconstrained, there may exist many 
low-mass, low- luminosity stars that can con- 
tribute significantly to the mass budget of 
a stellar population without creating a de- 
tectable increase in luminosity or change in 
color. For their derivation of the color-M/L 
relations, Bell & de Jong (2001) and Bell et 
al. (2003) adopted a truncated Salpeter IMF 
which derives from the constraint that bary- 
onic rotation curves should not over-predict 
observed rotation curves for spiral galaxies 
in Ursa Major (the Verheijen 1997 sample). 
With this constraint, they predict fewer low 
mass stars than a Salpeter IMF. These re- 
lations thus give an upper limit to the stel- 
lar mass present; a lower normalization of the 
color-M/L relation cannot be excluded. 

We re-derive this constraint on the upper 
limit to the IMF by re-calculating maximum 
disk fits for galaxies in the Verheijen (1997) 
sample without including a dark halo, as was 



6 




Fig. 3a. — Observed rotation curves (dots for H I, triangles for Ha or N II, thin dashed lines for models), 
baryonic rotation curves (solid lines), and dark matter rotation curves where applicable (thick dashed lines). 
For many of the Ha and N II rotation curves, the error bars arc smaller than the points at the resolution 
of the plots. The effects of the ±0.1 dex uncertainty in the color-M/L relations on the baryonic rotation 
curves are plotted as dotted lines. On the x-axis, the radii R25 and R = 2.2hm are marked with thick and 
thin bars, respectively. 
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Fig. 3b.— See Figure 3a. 
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Fig. 3c. — See Figure 3a. 
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Fig. 3d — See Figure 3a. 
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done in the original Verheijen (1997) fits. For 
these fits, we define maximal disk to be the 
greatest possible contribution of the bulge and 
disk to the observed rotation curve of the 
galaxy. To do this, we scale the LC-band de- 
rived stellar mass rotation curves as high as 
possible without over-predicting the observed 
rotation curves beyond the very inner parts. 
In order to best match the overall shapes of 
the observed rotation curves, we allow over- 
predictions of the inner 2 points of H I rota- 
tion curves due to beam-smearing and over- 
predictions of the inner parts of Ha rotation 
curves. We also re-scale the stellar and gas 
mass distributions from the distance used in 
Verheijen (1997) to 20.7 Mpc before fitting 
the maximal disks, unlike the approximate re- 
scaling done in Bell & de Jong (2001). Gas 
distributions derived from H I measurements 
are multiplied by 1.32 to account for helium, 
as in Verheijen (1997). In Figure 4, the result- 
ing maximum disk (M/L)* are plotted against 
the B — R colors of the galaxies corrected for 
extinction with the Tully et al. (1998) formal- 
ism. In this figure, we also plot maximum 
disk (M/L)* for those galaxies presented in 
this paper that have BRK photometry. All 
the points in Figure 4 are upper limits to the 
(M/L)*; galaxies cannot have values greater 
than those defined by the lower envelope in 
this plot without over-predicting their rota- 
tion curves beyond the very inner parts. It 
is apparent from Figure 4 that the Bell et al. 
(2003) color-M/L relation is consistent with 
all the galaxies to within ~ 0.1 dex, which is 
the uncertainty in the relations. 

Other than the IMF, there are five pri- 
mary sources of uncertainty in the determi- 
nation of the baryonic surface mass-density 
distributions: uncertainties in the distances 
to galaxies 2 , 0.1 dex spread in the color-M/L 




2 Of the 34 galaxies in our sample, 4 have distances 
estimated from Cepheid variables for which the un- 
certainties are much less than other types of distance 
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Fig. 4. — Maximum disk stellar (M/Lk)* versus 
reddening-corrected B — R colors for galaxies in 
the Verheijen (1997) sample (Xs) and in this paper 
(triangles). Plotted for reference are the color- 
M/L relations of Bell & de Jong 2001 (dashed 
line) and Bell et al. 2003 (solid line). 

relations, uncertainties in the zero-points of 
these relations, uncertainties in the photomet- 
ric zero-point calibrations, and uncertainties 
in the determination of a galaxy's inclination 
to the line of sight. Secondary sources of 
uncertainty generally have a small effect on 
the baryonic rotation curves. They include 
secondary photometric ( "bootstrap" ) calibra- 
tions, the neglect of interstellar gas content, 
and position angle uncertainties. Dust red- 
dening also plays a role in the uncertainty of 
the baryonic rotation curves; this will be dis- 
cussed in §4.3. Note that the transformation 
of Two Micron All Sky Survey (2MASS; Jar- 
rett et al. 2000; Cutrie et al. 2000; Jarrett et 
al. 2003) and SDSS photometry to the Kron- 
Cousins system introduces negligible uncer- 
tainties, as shown in Paper I. 

As an illustration, Figure 5a-d shows the 



measurements, and are typically ~ 10%. 
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effects of three of the primary sources of un- 
certainty on the baryonic rotation curve of 
NGC 157. In all the plots in Figure 5, we plot 
\Avb\/vf, where Avb is the difference between 
the velocity derived from the color-M/L rela- 
tions and this velocity affected by the named 
uncertainty. Figure 5a shows the effect of 
an uncertainty of ± 20% in the distance. 
The distance uncertainty causes a maximum 
change in \Avb\/vb of 0.24 which corresponds 
to a change in the velocity of nearly 40 km 
s _1 . The average change is \Avb\/vb = 0.09 
(|Ai>fc| ^ 13 km s _1 ). Figure 5b shows the 
effect of a systematic change of ±0.1 dex in 
the color-M/L relations. This uncertainty 
causes a maximum change in \Avb\/vb of 0.12 
(|A«5| ^25 km s^ 1 ), and an average change 
of 0.12 {\Av b \ ^ 20 km s" 1 ). In Figure 5c,d, 
the effects of introducing systematic errors 
in the photometric zero-point calibrations are 
shown. In Figure 5c, the effect of the actual 
zero-point uncertainties for NGC 157 is shown 
(a B = 0.03, a R = 0.03, o K = 0.04). This 
causes an average change in \Avb\/vb of 0.10 
(|Avb| ±2 17 km s -1 ). While most of our opti- 
cal photometric zero-points are good to < 5%, 
the worst uncertainty in an optical zero-point 
calibration for any of the surface brightness 
profiles that we use in this analysis is ±15%, 
as given in Table 3 of Paper I. All of the uncer- 
tainties on the near-infrared zero-points are 
~ 4%. To show the effect of a photomet- 
ric calibration that is not as good as that of 
NGC 157, in Figure 5d we show what would 
happen to the rotation curve if photometric 
zero-point errors were ±15% in the optical 
and ±4% in the near-infrared. Such uncer- 
tainties cause an average change in |Avb|/f;, 
of 0.13 (|Au 6 | ±2 22 km s" 1 ). 

In Figure 5e,f, we show the effects on the 
baryonic rotation curve of NGC 157 of chang- 
ing the inclination and position angle in the 
derivation of its surface brightness profiles. 
Typical uncertainties in position angles and 
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Fig. 5. — Effects of uncertainties on the baryonic 
rotation curve of NGC 157. For one side of the 
uncertainty (e.g., +0.1 dex), the |Avfc|/i>{, curve 
is plotted as a dotted line, and for the other side 
(e.g., -0.1 dex) the curve is plotted as a solid line. 
The effects of the following uncertainties are plot- 
ted: (a) ±20% uncertainty in the distance to the 
galaxy, (b) ±0.1 dex scatter in the color-M/L rela- 
tion, (c) actual photometric zero-point uncertain- 
ties for NGC 157, (d) photometric zero-point un- 
certainties for the galaxy with the poorest zero- 
point, (c) a change in the inclination of 5°, and 
(f) a change in the position angle of 5°. The radii 
2.2Hk and L25 are marked as thin and thick lines, 
respectively, on the x-axis in panel a. 

inclinations are both ~ ±5° (Paper I). In Fig- 
ure 5e, we plot the effect of a change in the po- 
sition angle of ±5°, and in Figure 5f we show 
the effect of a change in the inclination of ±5°. 
Uncertainties in the inclination have a larger 
effect (on average | ^2 0.03; \Avb\ — 5 km 
s -1 ) than those in the position angle (on av- 
erage \Avb\ — 0.007; \Avb\ — 1 km s -1 ). How- 
ever, both these uncertainties are small when 
compared to other sources of error discussed 
in this section. 

In Figure 6, to compare the color-M/L re- 
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baryonic rotation curves. 




r (kpc) 

Fig. 6. — Baryonic rotation curves derived from 
the color-M/L relations (plotted as dotted lines 
for the ±0.1 dex scatter in the relations) are com- 
pared with baryonic rotation curves derived from 
constant (M/Lk)* of 0.75 (thin solid line) and 1.0 
(thick solid line). The thin dashed line plotted for 
NGC 7606 is the model observed rotation curve of 
Courteau 1997. All other features of the plot are 
the same as Figure 3. 

lations to constant (M/Lk)*, baryonic rota- 
tion curves created with the color-M/L rela- 
tions and with (M/Lk)* of 0.75 and 1.0 are 
plotted for 4 example galaxies. The rotation 
curves created with (M/Lk)* of 0.75 and 1.0 
are approximately consistent with the ±0.1 
dex uncertainty of the color-M/L relations. 
This is similarly the case for all galaxies in 
our sample, and is due to the shallow slopes 
of the color-M/L relations (~ 0.14-0.18 in 
logioM/L). 

To summarize, other than the IMF, of the 
other sources of uncertainty in the determina- 
tion of the baryonic rotation curves, only the 
distance uncertainty, the ±0.1 dex scatter in 
the color-M/L relations, and zero-point un- 
certainties can produce non-negligible effects. 
These three sources of uncertainty vary from 
galaxy to galaxy and introduce scatter in the 



3.4. Effects Due to Dust 

Since the reddening vector in Figure 1 lies 
nearly parallel to the color-M/L relation, to 
first order, errors in foreground dust redden- 
ing estimates should not strongly affect the fi- 
nal relative derived masses of the stellar popu- 
lations, as foreground dust will systematically 
both redden and extinguish galactic light. In 
this section, we discuss the possible effects of 
dust reddening and extinction on absolute de- 
rived stellar masses. 

We examine the empirical inclination- 
dependent extinction correction of Tully et 
al. (1998), which should describe to first or- 
der the dust content of galaxies. We correct to 
face-on the total integrated colors and magni- 
tudes of galaxies in our sample. In doing this, 
we ignore radial dust gradients in the disk 
and the fact that dust reddening is likely to 
be large for the inner parts of galaxies. While 
the Lf-band magnitudes barely change when 
extinction corrected (0.05 mag on average), 
the B — R colors do (0.09 mag on average). 
This, however, should not affect the final de- 
rived stellar masses, since the color-M/L re- 
lations that we use have a very shallow slope 
in B — R color. To examine this effect further, 
for 4 galaxies that span the range of inclina- 
tions in our sample, we apply the Tully et al. 
(1998) extinction correction to their radial 
LT-band profiles, and use their extinction- 
corrected B — R colors to derive (M/Lk)*- 
From the resulting surface mass-density pro- 
files, we derive baryonic mass rotation curves. 
In Figure 7, we plot for the 4 galaxies the 
difference between the extinction-corrected 
baryonic rotation curve and the uncorrected 
one. For NGC 157, NGC 3726, NGC 7217, 
and NGC 7606, the rotation curves on aver- 
age differ by 1.3, 1.1, 2.7, and 5.0 km s _1 , 
respectively. The greatest difference is for 
NGC 7606 which has the largest inclination 
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Fig. 7. — Difference between the extinction- 
corrected (V) and uncorrected (V a ) baryonic ro- 
tation curves. The extinction-corrected baryonic 
rotation curves are derived from colors and magni- 
tudes corrected for extinction with the formalism 
of Tully ct al. 1998. 

of the four galaxies, 63.9°. These differences 
in velocity are all less than the uncertainties 
of the baryonic rotation curves themselves, as 
discussed in §3.3. 

Galaxy disks have been observed to have 
a relatively constant face-on dust opacity of 
~ 0.5 magnitude in the /-band (Holwerda 
et al. 2005). Correcting for this, Aj would 
increase the baryonic velocity by ~ 20 km 
s _1 . However, since the maximal disk limits 
in Figure 4 would also need to be corrected, 
and the color-M/L relations re-scaled accord- 
ingly, the net effect would more or less can- 
cel out. However, it is very likely that more 
complex dust models where the effects of dust 
depend on star/dust geometry are necessary, 
such as those by Disney (1989) and Gordon 
et al. (2001). We choose to leave this ap- 
proach to a future analysis, since the detailed 
radiative transfer codes are not publicly avail- 
able and are difficult to apply to real galaxies. 
Moreover, different galaxies may have vari- 
ous amounts of dust and different star /dust 
geometries which would cause dust to affect 



the baryonic rotation curve of each galaxy 
uniquely. 

4. Dark Matter 

For each galaxy, we compare its baryonic 
rotation curve derived in §3.2 with its ob- 
served rotation curve taken from the literature 
to derive a "dark matter rotation curve." 

4.1. Observed Rotation Curves 

Observed rotation curves are plotted in 
Figure 3; Table 1 lists the tracer and literature 
reference for each. Some galaxies have two ob- 
served rotation curves, one from Ha or N II 
observations which traces the inner parts, and 
another from H I observations which traces 
the outer parts. If errors for the observed ro- 
tation curves were given in the original ta- 
ble or plot that they were taken from, then 
these are used in this paper. If no errors are 
given, then we estimate them to be the dif- 
ference in rotation velocity between the ap- 
proaching and receding sides of the galaxy; ro- 
tation curves for which we estimate the errors 
in this manner are noted in Table 1. This esti- 
mate will generally give errors larger than the 
true measurement errors since it will be more 
affected by non-axisymmetric features such as 
spiral arms and slight warps in the gas distri- 
butions. Note that the 2-3 innermost points 
of the H I rotation curves and the outer few 
points of the Ha and N II rotation curves have 
a greater uncertainty than other points due 
to beam smearing and low signal-to-noise, re- 
spectively. For many of the rotation curves we 
have obtained data from the authors, but for 
a handful we could not. For those few galax- 
ies, rotation curves are extracted from plots 
in the literature with the DataThief program 
(Tummers 2000); these galaxies are noted in 
Table 1. Errors inherent to the extraction of 
a rotation curve vary from plot to plot, but 
tend to be < 5 km s _1 . 
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Galaxies marked with the reference "Math- 
ewson & Courteau" in Table 1 are rotation 
curves that were originally presented in Math- 
ewson, Ford, & Buchhorn (1992) and were 
later modeled by Courteau (1997). For these 
galaxies, we plot both the actual and model 
rotation curves in Figure 3. Since we are 
primarily interested in large-scale trends, we 
adopt the model rotation curves in the follow- 
ing analysis in order to avoid much of the fine 
structure inherent to the actual data. 

4.2. Dark Matter Rotation Curves 

At those radii where the observed rotation 
speed of a galaxy is greater than that of its 
baryons, the additional gravitational compo- 
nent is assumed to be due to dark matter. A 
dark matter rotation curve is derived as the 
square root of the difference of the squares of 
the observed rotation velocity and the bary- 
onic rotation curve velocity at each radius 
(Binney & Tremaine 1987). In doing this, it is 
assumed that the halos of galaxies are axially 
symmetric, the disk and halo are aligned, and 
the observed gas is in circular orbits. Dark 
matter rotation curves for the galaxies are 
plotted in Figure 3. 

For 10 galaxies, the baryonic rotation 
curves over-predict the observed rotation 
curves for some range in radius 3 . These 
galaxies are NGC 157, NGC 1559, NGC 2139, 
NGC2841, NGC 3198, NGC 4138, NGC 4698, 
NGC 5371, NGC 6300, and NGC 7083. The 
baryonic rotation curve for the -0.1 dex scatter 
in the color-M/L relations also over-predicts 
for the first six galaxies listed above. However, 
for these six galaxies, the over-prediction oc- 

3 We do not include in this count baryonic rotation 
curves that over-predict the inner parts of H I rota- 
tion curves, since they are affected by beam-smearing 
in this region. Also, Figure 4 suggests that this num- 
ber is less than f0, but we ignore the very inner parts 
of the galaxies when calculating maximal (M/Lk)*s 
in §3.1. 



curs only in the very inner parts where the ob- 
served rotation curve is affected by structures 
such as rings, bars, inner windings of spiral 
arms, and/or irregular morphology. Baryonic 
rotation curves are derived under the assump- 
tion of circular motion and need not trace such 
structures. For NGC 6300 and NGC 7083, 
the baryonic rotation curves for -0.1 dex do 
not over-predict the observed rotation curves. 
The galaxy NGC 6300 has a bar and a ring 
in the region where the over-prediction oc- 
curs, but the image of NGC 7083 shows no 
sign that its baryonic rotation curve should 
over-predict. This is somewhat acceptable, 
though, since the baryonic rotation curve for 
NGC 7083 for -0.1 dex does not over-predict. 

Two galaxies, NGC 4698 and NGC 5371, 
have baryonic rotation curves that over- 
predict their observed rotation curves, even 
for the -0.1 dex scatter in the color-M/L re- 
lations. For both these galaxies, the over- 
prediction cannot only be explained by bars 
and/or rings in the galaxy images. To ex- 
amine things further, we create baryonic ro- 
tation curves using the position angle and 
inclination that were used to derive their ob- 
served rotation curves; the resulting curves 
also over- predict. Next, we create baryonic 
rotation curves for the -20% uncertainty in 
the galaxies' distances and both the distance 
and color-M/L uncertainties; they are plot- 
ted in Figure 8. For NGC 4698, the baryonic 
rotation curve created by taking into account 
both these effects still over-predicts. However, 
for this galaxy, its distance was calculated 
with a Virgocentric infall calculation and was 
found to be triple- valued (Paper I). We derive 
a baryonic rotation curve for it with the dis- 
tance solution which is closer than the one 
chosen in Paper I (9.7 Mpc, as opposed to the 
chosen solution of 19.1 Mpc). This baryonic 
rotation curve is plotted in Figure 8 and over- 
predicts only the inner portion (~ 1.5 kpc) of 
the observed rotation curve. For NGC 5371, 
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Fig. 8. — Effects on the baryonic rotation curves 
of the -0.1 dex scatter in the color-M/L relations 
and the -20% uncertainty in distance for galax- 
ies where the baryons over-predict the observed 
rotation curve. Original baryonic rotation curves 
(thick solid lines), those derived by taking into ac- 
count the distance uncertainty (dashed lines), and 
those derived by taking into account both the dis- 
tance uncertainty and the color-M/L scatter (dot- 
ted lines) are triangles for Ha and small points for 
H I. For NGC 4698, we also plot the baryonic rota- 
tion curve derived from the adoption of a different 
distance to the galaxy (thin solid line). Observed 
rotation curves are plotted as points. The radius 
2.2hiR is marked as a thick solid line on the x-axis. 



the baryonic rotation curve created by taking 
into account the distance uncertainty over- 
predicts, while that created by taking into 
account both distance and color-M/L uncer- 
tainties does not. The sort of over-prediction 
that is observed for NGC 4698 and NGC 5371 
can also be due to such factors as a mis- 
measurement of the observed rotation curves 
and/or the effects of dust. A stellar popula- 
tion affected by dust may appear redder (and 
hence heavier and have a faster rotation) than 
it is intrinsically since the reddening effect of 
dust is slightly greater than its extinction ef- 
fect. These galaxies could also signal a need 
to lower the normalization of the color-M/L 
relations. If we did this, then many of the 
other galaxies in the sample would be sub- 
maximal, but the qualitative results of this 
paper would not change. 



4.3. Uncertainties in Dark Matter Ro- 
tation Curves 

Uncertainties in the dark matter rotation 
curves arise from a number of effects. The 
most significant are those inherent to the de- 
termination of the baryonic mass component, 
as discussed in §3.3. Other uncertainties in- 
clude: non-circular motions that perturb the 
underlying potential (i.e., spiral arms, bars, 
substructure) , statistical errors from the mea- 
surement of velocities in radial bins, system- 
atic errors in measuring the velocity (i.e., 
beam smearing and slit position angle er- 
ror) , and uncertainties in the measurement of 
the dynamical centers of the galaxies. There 
also may be differences between the centers 
of galaxies determined from photometry and 
those determined from the observed rotation 
curves. However, center measurements are 
not expected to differ much since photometric 
centers are always chosen to be the brightest 
pixel in the nucleus which coincides with the 
dynamical center of most galaxies. 

5. Maximal Disks 

The radius R = 2.2h is where the rotation 
curve of a self-gravitating exponential disk 
reaches its peak (Freeman 1970). A commonly 
used definition of a maximal disk is given by 
Sackett (1997) where the galaxy disk provides 
85% ± 10% of the total rotational support of 
the galaxy at 2.2/i#. In the following anal- 
ysis, we define a galaxy to have a maximal 
disk if it has a baryonic mass (disk and bulge) 
contribution to the observed rotation curve of 
> 90% at 2.2h, or similarly, a dark matter 
contribution to the observed rotation curve 
of < 10% at 2.2h. We choose to adopt this 
definition of maximal disk over that of Sack- 
ett (1997) since we perform detailed bulge- 
disk decompositions, and are therefore able 
to model the combined bulge plus disk bary- 
onic rotation curves. Furthermore, we use h 
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measured at K, and when if-band imaging is 
unavailable, we use H. Since near-infrared 
bands trace most of the mass of the stel- 
lar populations, a near-infrared scale-length 
is more analogous to the stellar mass scale- 
length of a disk. Disk scale-lengths measured 
at K are typically ~ 1.2 times shorter than 
those measured at B (de Jong 1996). 

In Figure 3, there are 4 galaxies that have 
submaximal disks, even if the -0.1 dex scatter 
in the color-M/L relations is taken into ac- 
count. Since the normalization of the color- 
M/L relations is an upper limit, the num- 
ber of galaxies in our sample that do not 
have maximal disks can in principle be much 
greater than 4. However, for those galaxies 
that we observe to have maximal disks, the 
overall shapes 4 of the inner parts of their ob- 
served rotation curves (within ~ R25) are gen- 
erally matched by the shapes of their baryonic 
rotation curves. This is evidence that many 
of the disks have a significant baryonic com- 
ponent in their inner parts. 

These 4 galaxies are NGC 3319 5 , NGC 3992, 
NGC4062, and NGC 7606; they have bary- 
onic mass contributions to their total masses 
at 2.2h IR of 22%, 50%, 57%, and 68%, re- 
spectively. In Figure 9 we plot for these 4 
galaxies the baryonic rotation curves with 
and without uncertainties in distances and the 
color-M/L relations taken into account; these 
baryonic rotation curves still under-predict 
the observed rotation curves for NGC 3319, 
NGC 3992, and NGC 4062. However, for 

4 This is distinguished from the "bumps and wiggles" 
of the observed rotation curves that are likely due to 
small-scale features such as bars and spiral arm per- 
turbations (Palunas & Williams 2000; Kranz, Slyz, & 
Rix 2003; Slyz, Kranz, & Rix 2003). 

5 For NGC 3319, we do not take into account its in- 
nermost 2 observed rotation curve points since they 
are measured from H I and are likely affected by 
beam-smearing. Also, the distance measurement for 
NGC 3319 is from Cepheid variable stars and has an 
uncertainty of only ~ 10%. 



NGC 7606, the baryonic rotation curve that 
takes into account both these effects follows 
the observed rotation curve fairly well. For 
NGC 4062, its distance was calculated us- 
ing a Virgocentric infall calculation and was 
found to be triple- valued (Paper I). We cre- 
ate baryonic rotation curves for NGC 4062 
with the other 2 distance solutions (17.6 and 
24.4 Mpc), both of which are greater than 
the one chosen in Paper I; they are plotted 
in Figure 9. These baryonic rotation curves 
have a strange behavior since they under- 
predict the observed rotation curve in the 
inner parts, but match it in the outer regions. 
This may be the effect of such factors as an 
underestimate of the stellar mass in the inner 
parts of the galaxy due to blue spiral arms 
in the 5-band image (see image in Figure 1 
of Paper I), a poor bulge/disk decomposition, 
and/or an underestimate of the uncertainties 
of the observed rotation curve. In summary, 
while NGC 4062 and NGC 7606 may not have 
submaximal disks, NGC 3319 and NGC 3992 
likely do. 

There are a few galaxies that have marginally 
submaximal disks. For these galaxies, the 
baryonic rotation curves for the +0.1 dex 
scatter in the color-M/L relations result in 
maximal disks. These galaxies are NGC 1241, 
NGC 2139, and NGC 2280, and they have 
baryonic mass contributions to their total ro- 
tational support at 2.2h[R of 75%, 73%, and 
72%, respectively. 

One striking example of a galaxy that is 
likely close to maximum disk is NGC 157. 
This galaxy has Ha and H I rotation curves 
that have a sudden steep decline at ~ 8kpc 
(~ 3hx) and flatten afterward, which results 
in a hump-like structure. While this peculiar 
behavior cannot be absolutely confirmed by 
Ryder et al. (1998), there are strong lines of 
evidence presented in their paper that point 
to this hump-like structure as physical. This 
same structure is also found in the baryonic 



17 







| "" | " i ' | ' i ' i | i i i i | i ' i ' | i- 




400 




i ' i | i ' i ' | ' i ' i | i ' ' ' | i '. 




150 


- NGC3319 1 








NGC 3992 






: - Iir ***j 


CO 


300 




t i 1 1 1 




100 










" i i i i ■ 


s 






6 


200 






















50 






100 






> 






> 










n 


i , , , , 1 — i , , , , i , , , , i — i ," 








lllllllllMlllKlll 



5 10 15 20 25 10 20 30 40 




02466 5 10 15 



r (kpc) r (kpc) 

Fig. 9. — Same as Figure 8, but for galaxies where 
the baryons under-predict the observed rotation 
curves, and for the +0.1 dcx scatter in the color- 
M/L relations and the +20% distance uncertainty. 
For NGC 4062, baryonic rotation curves derived 
for 2 different distance determinations are plotted 
as thin solid lines. 

rotation curve for NGC 157. Moreover, a dark 
matter halo in the shape of a NFW model 
is not consistent with this structure unless 
baryons make a significant contribution to the 
inner ~ 5 scale-lengths of NGC 157. In sum- 
mary, although a lower normalization of the 
IMF cannot be ruled out, it seems unlikely 
that at least NGC 157 is strongly submaximal. 

6. Dark and Baryonic Matter Scaling 
Relations 

In this section, we examine scaling relations 
for dark and baryonic matter. We use quanti- 
ties from Tables 1 and 4 of Paper I, and derive 
others: V tot . max the maximum observed rota- 
tion curve velocity, Vb,max the maximum bary- 
onic rotation curve velocity, R(Vb,max) the ra- 
dius at which V = Vb,max-, and Mb the bary- 
onic mass. These quantities are listed in Table 
2, and if they are derived from the baryonic 



mass distributions, we also tabulate the dif- 
ferences in their values for renormalizations 
of the color-M/L relations by ±0.1 and -0.3 
dex. We choose to create the quantity Vb,max 
because it can be derived from imaging alone 
and obviates the need for much more expen- 
sive and time-consuming line width observa- 
tions needed to obtain Vt t,max- For bright 
spirals, Vj, ,max should not differ from Vtot,max 
by very much, and due to the flat and usually 
noisier nature of the observed rotation curves, 
the measurement of Vb,max is more straightfor- 
ward than that of Vtot,max ■ This quantity will 
be discussed further in §6.2. 

We derive two quantities from the dark 
matter rotation curves calculated in §4.2: i? 10 , 
the radius where dark matter contributes 10% 
to the velocity of the observed rotation curve, 
and Rx, the radius where the dark matter 
contribution equals that of the baryons (the 
"cross-over radius" ) . The quantity i?io is sim- 
ilar to i?mD of Salucci (2001) and R t of Persic 
et al. (1996). The radius Rx is analogous to 
R 2 -.i of McGaugh & de Blok (1998). The radii 
Rio and Rx are listed in Table 2 along with 
the difference in values for renormalizations 
of the color-M/L relations by ±0.1 and -0.3 
dex. Note that for some galaxies we cannot 
measure i?io or Rx ; this is in general because 
their rotation curves do not extend far enough 
in radius. 

6.1. Baryonic Scaling Relations 

In Figure 10, basic physical parameters of 
galaxies are plotted versus M&: total abso- 
lute B and K-b&nd magnitudes, Hubble T- 
type, R25, and B and iv-band central surface 
brightnesses (fJ. ,B, ^o,k)- They have the fol- 
lowing correlation coefficients: 0.91 (B), 1.00 
(K), 0.40 (Hubble T-type), 0.78 {R 25 ), 0.29 
(ij- ,b), and 0.58 {^ ,k)- Integrated magni- 
tudes, sizes, fi ,Ki and maximum rotation ve- 
locities correlate very well with Mb, such that 
galaxies with greater Mb are brighter, larger, 
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Fig. 10. — Basic physical parameters of galax- 
ies versus baryonic mass. Open circles represent 
galaxies from the SDSS that have partial imaging. 
A typical error bar for Mb is plotted in panel a; 
error bars for the other quantities are the size of 
the data points at the resolution of the figure. To 
make the plot in panel c such that points do not 
overlap, we add fractions with values < 1 to the 
integer Hubble T-types. 

and have brighter fJ, ,K- Morphology and (i 0> b 
also correlate with M b , but to a lesser extent. 

In Figure 11, relations are plotted for 

.max 

and logioVtot,max, which trace 
the baryonic and total mass components of 
the galactic potential, respectively. In Fig- 
ure 11a, along the lines of Roberts & Haynes 
(1994), we plot log w V to t,max versus log w L B 
and find good agreement with their results. 
In Figure lib, we plot log 10 V b:Tnax versus 
logioVW, max] they are strongly correlated 
with a correlation coefficient of 0.82. These 
velocities should be equivalent for maximal 
disks, which many of the galaxies in our sam- 
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Fig. 11. — Relations for V tot , m ax and V btTnax , 
which trace the total and baryonic mass compo- 
nent of the galactic potential, respectively. Open 
circles represent galaxies from the SDSS that have 
partial imaging. Error bars for T-type, Lb, and 
i?25 are the size of the data points at the resolution 
of the figure. 

pie have. In addition, a tight correlation indi- 
cates that there is not a wide spread in the de- 
gree of maximality. The quantity logio Vtot, max 
also correlates with Hubble T-type and R25 
such that galaxies that rotate faster have ear- 
lier T-types and are larger. The quantity 
R(V b , max) correlates with logio-^6 with a cor- 
relation coefficient of 0.66, but correlates rel- 
atively weakly with logioV& mox , logi V to t ,max, 
and Ho^k with correlation coefficients of 0.08, 
0.35, and 0.16, respectively. 

As discussed above, the maximum rotation 
speeds predicted from the baryon distribu- 
tions are tightly correlated with the observed 
maximum rotation speeds. To examine this 
relation further, we plot the ratio of V\,^ max to 
V to t,max versus [i 0)K in Figure 12, and a find 
root mean square (rms) deviation from unity 
of only 0.18. Using this result, one can cre- 
ate a Tully-Fisher relation from two passband 
surface photometry and a redshift alone (e.g., 
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Fig. 12. — Ratio of maximum baryonic veloc- 
ity to maximum observed velocity versus H ,k- 
SDSS galaxies with partial imaging are plotted as 
open circles. Outliers are labeled: NGC 4062 and 
NGC3992 have submaximal disks, NGC 4138 has 
evidence of kinematic disturbance, and NGC 3319 
has both a submaximal disk and a kinematic dis- 
turbance. 

by using the SDSS and the color-M/L rela- 
tions for g — r and (M/L r ) m ). The 4 outliers 
in Figure 12 are galaxies that have submaxi- 
mal disks and/or are kinematically disturbed. 
The relation between Vb,max and Vtot,max is 
not expected to be as tight for less luminous 
spirals, but could possibly be calibrated for 
such a population. 

In Figure 13a, b we plot the baryonic 
Tully-Fisher relations for V to t,max and V btmax . 
Weighted bisector fits are given in the form: 
Mb = M2ooV a where M200 is the stellar mass 
in units of 1O 1O M0 for a galaxy on the Tully- 
Fisher relation with V = 200 km s _1 . They 
are: M 200 = 5.40 ± 0.03 and a = 3.4 ± 0.3 
with a rms residual of 0.20 for Vt t,max, and 
M200 = 6.98 ± 0.03 and a = 2.8 ± 0.3 with 
a rms residual of 0.22 for Vb,max- If we 



Fig. 13. — Baryonic Tully-Fisher relations for 
y tot. max and V b , max , and the sizc-baryonic mass 
relation are plotted in panels a-c. Weighted bi- 
sector fits with and without NGC 3319 are plotted 
as black and gray lines, respectively. In panel d, 
the V to t,max residual of the Tully-Fisher relation 
is plotted versus the size residual from the size- 
baryonic mass relation. Galaxies from the SDSS 
that have partial imaging are plotted as open cir- 
cles. 

do not include the outlier in these two re- 
lations, NGC 3319, we obtain the following 
fits: M 200 = 5.36 ± 0.03 and a = 3.3 ± 0.3 
with a rms residual of 0.19 for Vt t,max, and 
M200 = 6.75 ± 0.03 and a = 3.1 ± 0.3 with 
a rms residual of 0.21 for Vb jma x- The galaxy 
NGC 3319 is both kinematically disturbed on 
one side (Moore & Gottesman 1998), and has 
the smallest stellar mass in our sample. For 
most of the galaxies in our sample, we do 
not have H I or molecular hydrogen measure- 
ments. If gas masses were included, the Tully- 
Fisher slope of Mb on Vtot,max would flatten, 
as lower mass galaxies (e.g., NGC 3319) have 
a larger gas content than more massive galax- 
ies (e.g., Verheijen 1997). 
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We refrain from comparing the slopes 
found here and in the literature (e.g., Bell & 
de Jong 2001; McGaugh 2005; Pizagno et al. 
2005) for the following reasons. Such fits are 
sensitive, among other things, to the velocity 
measurements used (i.e., Vtot,maxi Vb.maxt Vf lat ,V 2 .2, 
W20), the range in stellar mass of the sam- 
ples, and the galaxy samples used (range 
in Hubble type, field versus cluster environ- 
ments). In particular, Gurovich et al. (2005) 
find a break in the stellar mass and baryonic 
Tully-Fisher relations for less massive galaxies 
(10 7 - 6 < M* < 1O 9 - 6 M ). This break is such 
that the Tully-Fisher relation for less mas- 
sive galaxies is steeper than that for galaxies 
with stellar masses from ~ 1O 9,6 -1O 1L2 M . 
Gurovich et al. (2005) find slopes for the lower 
and upper mass ranges of their fZ-band sam- 
ple, using W20 as a velocity measurement, 
of 4.4 ± 0.3 and 3.3 ± 0.3, respectively. This 
break/curvature makes a comparison between 
slopes of different samples risky at best, es- 
pecially when different selection criteria are 
used. A 2D Kolomogrov-Smirnov test would 
be ideal, but beyond the scope of this work 
considering the modest stellar mass range of 
our sample. 

The size-baryonic mass relation, where size 
is parameterized by hjR, is plotted in Fig- 
ure 13c. We perform a weighted bisector fit 
and find h IR = 0.67 ± 2.31 (M 6 /lO lo M ) a6O±a08 
kpc with a rms of 0.16. We compare the 
Vtot,max residuals of the Tully-Fisher relation 
with the size residuals of the size-mass rela- 
tion in Figure 13d, and find a scatter plot. If 
galaxy disks are maximal and the dark matter 
halos of galaxies of different surface brightness 
are identical, then it is predicted by McGaugh 
& de Blok (1998) and Courteau & Rix (1999) 
that there should be a correlation between 
these residuals. Courteau & Rix (1999) ar- 
gue that the lack of residual correlation, es- 
sentially the surface brightness independence 
of the Tully-Fisher relation, implies that all 



disks are submaximal. This argument, how- 
ever, is not entirely straightforward, since in 
the Pizagno et al. (2005) model for disk col- 
lapse within a halo, variations in dark matter 
halo parameters can create enough scatter 
in the Tully-Fisher relation to hide the pre- 
dicted correlation. Furthermore, McGaugh 
& de Blok (1998) and Sellwood (1999) ar- 
gue that submaximal disks do not solve the 
problem, and one is left with a fine-tuning 
problem: either disk M/Ls correlate with 
surface brightness, halo contributions vary, 
or Newtonian dynamics falters. Figure 12 il- 
lustrates this fine-tuning problem by showing 
that V btmax /V to t,max is independent of surface 
brightness, so that the total velocity "knows 
about" the baryonic contribution. Changing 
the normalization of the color-M/L relations 
re-scales the y-axis of the plot, but the ratio 
remains independent of surface brightness. 

6.2. Dark Matter Scaling Relations 

In Figures 14a-d and 15, the baryon/dark 
matter equality radius in units of the disk 
scale length, Rx/hm, is plotted versus a num- 
ber of galactic parameters: V to t,max, V bjmax , 
Hubble T-type, M b , B, K, fi ,B, and fJ> ,K- In 
Figure 14e, we plot Rx versus hiR. There are 
fewer points in these figures than in previous 
ones since not all galaxies have a dark mat- 
ter rotation curve that allows us to determine 
Rx- The correlation coefficients for these re- 
lations and for R(Vb y max) with Rx are listed in 
Table 4. The radius Rx is found to correlate 
most strongly with V b) max and very strongly 
with M b , Hubble T-type, and V to t,max- There 
are no changes in the relative strengths of the 
relations if the color-M/ L relations are renor- 
malized by ±0.1 or —0.3 dex. 

The trends with Rx that we find are qual- 
itatively in agreement with those in the lit- 
erature (e.g., de Blok & McGaugh 1997; Mc- 
Gaugh k de Blok 1998; Pizagno et al. 2005), 
with some exceptions that are discussed here. 
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Fig. 14. — Relations for Rx/hm and Rx- In panels a-e, triangles represent galaxies from SDSS that 
only have partial imaging. Error bars represent the values of the quantities derived by renormalizing the 
color-M/L relations by ±0.1 dcx; upper and lower limits to these values are plotted in gray. To make the 
points in panel c such that they do not overlap, we add fractions with values < 1 to the integer Hubble 
types. In panel f, galaxies are plotted with different symbols according to their hm: those with hm < 1 kpc 
as open circles, those with 1 kpc < hiR < 3 kpc as filled circles, and those with 3 kpc < Iiir as open 
triangles. The galaxies NGC 2841, NGC 3319, and NGC 4138 are labeled in all panels; they have evidence 
of kinematical disturbances despite their normal appearances. The three outliers in Figures 14 and 15 
are NGC 2841, NGC 3319, and NGC 4138. Despite their normal optical morphologies, these galaxies have 
signatures of kinematic disturbances. The galaxy NGC 2841 has a warp in its outer H I disk (Bosma 1981) 
and an indication of a counter-rotating stellar component for 5" < r < 12" (Silchenko, Vlasyuk, & Burenkov 
1997). Similarly, NGC 4138 has both a counter-rotating disk and a significant warp in its outer H I disk 
(Jore, Broeils, & Haynes 1996). This warp may be the cause of the decline of its H I rotation curve, and 
hence what makes this galaxy an outlier. The third outlier, NGC 3319, is discussed in §6.1. 
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Fig. 15. — Relations for Rx/hm and inte- 
grated magnitudes and central surface bright- 
nesses. Symbols and error bars are the same as 
in Figure 14. The galaxy NGC 3319 has a fj, 0iK 
measurement, but not a K measurement. Error 
bars along the x-axis are approximately the point 
size at the resolution of the plot. 

McGaugh k de Blok (1998) and de Blok k 
McGaugh (1997) find that fi 0: B plays a major 
role in relations with galaxy properties, in- 
cluding with i?2:i, which is analogous to our 
Rx parameter. We find that Rx/hm is cor- 
related with fi ,B, but that it is even more 
correlated with [x ,K- We do not find a corre- 
lation between the mass discrepancy and // ,K 
(in Figure 12 we plot the inverse mass discrep- 
ancy), but McGaugh k de Blok (1998) do find 
a relation between total M/Lb at Ah and \x ,B- 
Also, we find that hiR correlates strongly with 
Rx-, once outliers are removed, which is not 
in qualitative agreement with other studies: 
McGaugh k de Blok (1998) and Zavala et al. 
(2003) find that h does not correlate with to- 
tal Mj L b evaluated at Ah or the mass discrep- 
ancy evaluated at maximum rotation velocity, 
respectively. However, the results of Pizagno 
et al. (2005) for relations with h are consistent 



with ours; they find that more compact galax- 
ies have a larger mass discrepancy (measured 
at 2.2/ij) than larger galaxies. Zavala et al. 
(2003) find a number of relations involving a 
quantity similar to the mass discrepancy that 
are in contradiction to those presented in this 
paper and in the literature. 

Here we discuss possible causes of these dis- 
crepancies. Our differences with McGaugh k 
de Blok (1998) may be traced to the quanti- 
ties used in the analyses: they use M/Lb and 
hs, while we use Rx and hiR. Luminosity in 
the .B-band should have more scatter in re- 
lations with galaxy properties than baryonic 
mass derived from a combination of optical 
and near-infrared data. This is because the 
-B-band is more affected by star formation and 
extinction, while the If -band is a better tracer 
of stellar mass. In addition, a near-infrared 
scale-length is more analogous to the bary- 
onic mass scale-length of a disk, and hence 
should have less scatter in its relations with 
galaxy properties than an optical scale-length. 
Adding some credence to the hypothesis that 
different quantities are the root of the discrep- 
ancy, McGaugh k de Blok (1998) find that 
relations between R2.1 and both \x 0) b and Mb 
are such that brighter galaxies have larger val- 
ues of R2.1, consistent with our results. That 
McGaugh k de Blok (1998) find a correlation 
between M/Lb evaluated at Ah and ii ,b is 
likely also due to the larger range in surface 
brightness in their sample. 

In Figure 14f, i?io is plotted versus Rx, and 
galaxies are plotted as different symbols ac- 
cording to their value of h[R. Although there 
is some scatter, as hjR increases, these two 
radii move further out in the disks in tandem. 
That there is such a tight relation between 
Rio and Rx tells us that dark matter contri- 
butions to the observed rotation curves must 
increase in a characteristic way between these 
two radii for all galaxies. This is likely due 
to the combined effects of quasi-exponential 
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disks and observed rotation curves that are 
nearly flat. If this is correct, then Rio should 
have more scatter in its relations with galaxy 
properties than Rx because the observed ro- 
tation curves should not yet be flat in the re- 
gion where R\o is measured, which is found to 
be the case. 

7. Radial Behavior of Dark Matter 

In Figure 16, the dimensionless parame- 
ter P(r) = M b (r)/M tot (r), where M b is the 
baryonic mass, and its inverse are plotted for 
galaxies with an appreciable dark matter con- 
tribution. It measures the fractional contri- 
bution of baryons to the gravitational poten- 
tial function of radius in a galaxy, and 
is akin to the /3 parameter defined by Salucci 
(2001) and similar parameters used in many 
other papers, except here it is evaluated at 
all radii. Where (3{r) = 1, the baryonic 
mass of a galaxy accounts for its observed 
rotation curve, j3(r) = 0.5 at r = Rx, and 
P(r) = where the dark matter accounts 
for its observed rotation curve. If the bary- 
onic model over-predicts the observed rota- 
tion curve, then (3{r) > 1. In Figure 16a, (3(r) 
is plotted versus radius in terms of hm. For 
many galaxies, the observed rotation curves 
are entirely accounted for by baryons in the 
inner parts. Beyond this region, baryonic 
mass falls off as dark matter begins to dom- 
inate. Other galaxies are dark matter domi- 
nated throughout. Curves in Figure 16a have 
different line types that correspond to ranges 
of Vb : max- For most of the galaxies, as Vb.max 
increases, so does the proportion of baryonic 
to dark matter at all radii such that the fastest 
rotators are observed to be dominated by 
baryons until quite far out into their disks. 
However, there is clearly much scatter about 
this trend; this conclusion can also be inferred 
from Figure 14f. 

In Figure 16b, (3(r) is plotted versus radius 
in units of Rx- This choice of radial coordi- 
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Fig. 16. — (3{r) = M b (r) / M to t{r) and its inverse 
for galaxies with an appreciable dark matter con- 
tribution. (3(r) is plotted versus r/hm and r/Rx 
in parts a and b, respectively; its inverse is plotted 
versus r/Rx in panel c. Different line types corre- 
spond to the Vb t max of the galaxies: Vb, ma x > 250 
km s" 1 (solid line), 201 < V b>max < 250 (dot- 
ted line), 120 < Vb, max < 201 (dashed line), and 
Vb, max < 120 (dot-dash line line). 

nate causes the curves to overlap at r = Rx 
(where (3(r) = 0.5), and thus allows for a bet- 
ter comparison of their radial behavior. In 
Figure 16c, we plot the somewhat less intu- 
itive function (3~ l {r) since it increases linearly 
with r /Rx ■ There is little variation in the be- 
havior of the curves in Figures 16b, c. Such 
regularity can be explained in terms of quasi- 
exponential disks and flat rotation curves as 
follows: The rotation curve of a galaxy is 
nearly flat beyond r ~ 2hjR, or it is at least 
a slow function of r and varies less than V&(f) 
in this region. In addition, much of the bary- 



24 



onic mass of a galaxy is enclosed at r ~ 2hm, 
which causes the baryonic rotation curves to 
be roughly Keplerian (14(r) oc 1/r) beyond 
this radius. Therefore, for all galaxies, beyond 
~ 2hm, it should be the case that (3{r) = 
V 2 /V 2 ot oc 1/r, which is what we observe. All 
the (3(r) curves overlap nicely when plotted 
versus radius normalized to Rx since this ra- 
dius is located beyond 2hm, and is generally 
in the falling part of the baryonic rotation 
curves. We parametrize the trend of with 
r/Rx in a simple universal manner that holds 
for all galaxies: /T 1 = 1.71(r/R x ) + 0.021. 

This relation is consistent with the predic- 
tion of Palunas & Williams (2000) that ei- 
ther the contribution of dark matter within 
the optical radius of galaxies is small or that 
the distribution of dark matter is coupled to 
that of the luminous matter. This relation 
is also similar in spirit to the parameteriza- 
tion of the mass discrepancy-acceleration re- 
lation of McGaugh (2004). Since (3~ l is the 
mass discrepancy, and a = V 2 /r oc 1/r for flat 
rotation curves, a natural mass discrepancy- 
acceleration relation arises. 

8. Comparison With Theories of Halo 
Density Distributions 

We compare the derived dark matter pro- 
files with an analytical function designed to 
parametrize the density profile of dark matter 
halos in N-body simulations. In particular, 
we compare our data to the NFW formulation 
for dark matter halo density profiles (Navarro, 
Frenk, & White 1996) with and without tak- 
ing into account adiabatic contraction of the 
halos. To do this, for each galaxy, we fit its 
observed rotation curve with a total mass ro- 
tation curve created from the addition of its 
baryonic rotation curve to a grid of halo mod- 
els, with a reduced x 2 statistic. Uncertain- 
ties in the observed rotation curves are taken 
to be 10 km s -1 ; results do not differ signifi- 
cantly if the error bars plotted in Figure 3 are 



used. To test the most common implementa- 
tion of dark matter contraction (Blumenthal 
et al. 1986), we adiabatically contract the grid 
of halo models according to the radial density 
distribution of baryons in the galaxies follow- 
ing the formalism of Dutton et al. (2005), and 
perform the fits again. The NFW fitting for- 
mula has 2 free parameters that we fit for 
using a grid covering: 0.5 < C200 < 20 and 
0.5 V tot ,max < V 200 < 2.5 V tot ,max ■ 

In Figure 17, we plot the best-fit halo mod- 
els for 8 example galaxies, and in Table 3 
we list the best-fit parameters along with the 
reduced x 2 of all the fits. Most of the fits 
are very poor; the average reduced \ 2 values 
with and without adiabatic contraction are: 
7.1 and 4.2 for the original color-M/L rela- 
tions, 14.0 and 7.3 for the +0.1 dex renormal- 
ization of the color-M/L relations, 4.8 and 3.4 
for —0.1 dex, and 3.9 and 4.3 for —0.3 dex. We 
do not perform fits for galaxies when most of 
the observed rotation curve is accounted for 
by the baryonic rotation curve. 

For baryonic rotation curves derived from 
the original color-M/L relations, nearly all fits 
to the NFW models without contraction have 
a smaller reduced x 2 than those where con- 
traction was performed. The exceptions are 
NGC4062 and NGC3992, which likely have 
submaximal disks. For galaxies that have a 
large baryonic contribution to the inner parts 
of their rotation curves, adiabatic contrac- 
tion over-contracts the inner parts of their 
dark matter halos such that the resulting to- 
tal mass rotation curves over-predict the mea- 
sured rotation curves, as in Weiner, Sellwood, 
& Williams (2001), for example. For fits to 
the NFW models where the baryonic rotation 
curves were derived from color-M/L relations 
that were renormalized by +0.1 dex, the situa- 
tion is exactly the same. For -0.1 dex, 7 galax- 
ies have better fits when contraction is taken 
into account, and for -0.3 dex, 15 galaxies have 
a better fit. This implementation of dark mat- 



25 




Fig. 17. — NFW halo fits for 8 example galaxies without (left) and with (right) adiabatic contraction. 
Best-fit NFW models are plotted as dashed lines, baryonic rotation curves as thick solid lines, total rotation 
curves (sums of the best-fit NFW models and the baryonic rotation curves) as thin solid lines, and observed 
rotation curves as points. For those galaxies with a rotation curve from Mathewson et al. 1992 that is 
modeled by Courteau 1997, we only plot the model as a dotted line for clarity. 
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Fig. 18. — Best-fit parameters for NFW fits with 
and without adiabatic contraction (filled and open 
symbols, respectively) . The median and the 68 per 
cent of the C200 values measured in numerical N- 
body simulations of dark matter halos are plotted 
as solid and dashed lines, respectively. When the 
best-fit halo is at the edge of the parameter space 
searched, it is plotted as a triangle, otherwise it 
is plotted as a circle. In part a, best-fit parame- 
ters for the original and ±0.1 dex renormalizations 
of the color-M/L relations are plotted; in part b, 
best-fit parameters for the -0.3 dex rcnormaliza- 
tion are plotted. 

ter contraction works best when baryons do 
not account for all of the inner parts of the 
observed rotation curve. 

In Figure 18, we plot the best-fit values 
of V200 and C200 and the range of these pa- 
rameters found in N-body simulations (Bul- 
lock et al. 2001; Eke, Navarro, & Steinmetz 
2001). Even under the assumption of no adia- 



batic contraction, the derived concentrations 
for these halos are low compared to those 
measured from simulations. The concentra- 
tions for —0.3 dex are in slightly better agree- 
ment with simulations, but are still quite low. 
In concert with this, McGaugh, Barker, h, 
de Blok (2003) and McGaugh (2004) found 
concentrations for NFW halos that are too 
low for a standard ACDM universe. In ad- 
dition, Alam, Bullock, k Weinberg (2002) 
examined low surface brightness galaxies and 
found their dark matter halos to be under- 
concentrated compared to what is expected 
from a standard ACDM universe, even if 
they assumed the galaxies to be dark matter- 
dominated. 

9. Summary 

We have decomposed the rotation curves 
of 34 nearby bright spiral galaxies into bary- 
onic and dark matter components by applying 
color-M/L relations to near- infrared and op- 
tical photometry, and find the following: 

• The dark-to-luminous matter distribu- 
tions are self-similar once scaled by 
Rx, the radius where the baryonic and 
dark matter contributions to the rota- 
tion curve are equivalent. This behav- 
ior is parameterized by a simple func- 
tion whose form is due to the quasi- 
exponential nature of galaxy disks and 
rotation curves that are nearly flat after 
an initial rise. This result is indepen- 
dent of the normalization of the color- 
M/L relations. 

• The radii i?io and Rx, where dark 
matter contributes 10% and 50% to 
the rotation of galaxies, respectively, 
correlate with galaxy properties. The 
strongest correlation with Rx is with 
the maximum baryonic rotation speed 
such that galaxies with Rx measure- 
ments that lie further out in their disks 
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rotate faster. The next strongest corre- 
lations are equivalently baryonic mass, 
observed rotation speed, and Hubble T- 
type. Contrary to what is expected from 
previous studies (de Block & McGaugh 
1998; Zavala et al. 2003), L?-band cen- 
tral surface brightness is not found to 
be the main driving force in relations 
with Rx for this sample of bright galax- 
ies. The radii Rx and R\$ move out in 
galaxy disks in tandem, consistent with 
a self-similar dark-to-luminous matter 
distribution among galaxies. 

• We confirm the normalization of the 
color-M/L relations given in Bell et al. 
(2003), which is analogous to an upper 
limit on the IMF based on the dynamics 
of disk galaxies. A more careful analy- 
sis of the data used in both Bell & de 
Jong (2001) and Bell et al. (2003) is per- 
formed, and data from this paper are 
included in the sample. 

• All but 4 of the 34 galaxies in our sample 
are close to maximal disk. Two of these 
4 galaxies have maximal disks within 
uncertainties. A prime example of a 
galaxy that cannot have a submaximal 
disk given the current formulations of 
dark matter halos is NGC 157 which has 
a pronounced hump-like structure in its 
observed and baryonic rotation curves. 

• Maximum rotation velocities predicted 
from the baryon distributions of galax- 
ies (Vb,max) are tightly correlated with 
the observed maximum rotation speeds 
(Vtot ,max)- Using this, a baryonic Tully- 
Fisher relation can be created based on 
two-passband surface photometry and 
a redshift alone (e.g., with the SDSS 
and color-M/L relations for g — r and 
(M/L r )*, or for more distant redshift 
surveys with Hubble Space Telescope 
imaging). Such a relation could possibly 



be calibrated for lower luminosity galax- 
ies. 

• We find generally poor fits for the NFW 
parameterization for dark matter ha- 
los due to the significant baryonic con- 
tributions found in the inner parts of 
most galaxies. The concentrations of 
the best-fit NFW halos are therefore 
much lower than what is expected for 
galaxies in a standard ACDM universe. 
This is even the case when the color- 
M/L relations are renormalized by -0.3 
dex. In order to have better fits, a 
normalization even lower than -0.3 dex, 
where baryons contribute very little to 
the total mass in the inner parts of ro- 
tation curves, would have to be imple- 
mented. Adiabatic contraction, as it is 
normally implemented, makes these fits 
worse by adding more dark matter to 
the inner parts of galaxies. 
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Table 1 



Galaxy 


Tracer 3, 


Reference(s) 


NGC 157 


FP Ha, H I 


Fridman ct al. 2001, Ryder et al. 1998 b 


NGC 289 


HI 


Walsh, Stavelcy-Smith, & Oosterloo 1997 b 


NGC 488 


Ha 


Peterson 1980 b < c 


NGC 908 


Ha 


Mathewson, Ford, & Buchhorn 1992 


NGC 1087 


Ha 


Rubin et al. 1985 


NGC 1090 


Ha, H I 


Courteau 1997, Gentile et al. 2004 


NGC 1241 


Ha 


Mathewson, Ford, & Buchhorn 1992 


NGC 1385 


Ha 


Mathewson, Ford, & Buchhorn 1992 


NGC 1559 


Ha 


Mathewson, Ford, & Buchhorn 1992 


NGC 1832 


Ha 


Mathewson, Ford, & Buchhorn 1992 


NGC 2090 


Ha 


Mathewson, Ford, & Buchhorn 1992 


NGC 2139 


Ha 


Mathewson, Ford, & Buchhorn 1992 


NGC 2280 


Ha 


Mathewson, Ford, & Buchhorn 1992 


NGC 2841 


FP Ha, H I 


Blais-Ouellette et al. 2004, Giraud 1998 c 


NGC 3198 


FP Ha, H I 


Corradi et al. 1991 b , van Albada et al. 1985 


NGC 3223 


Ha 


Mathewson, Ford, & Buchhorn 1992 


NGC 3319 


H I 


Moore & Gottesman 1998 b ' c 


NGC 3521 


H I 


Sanders 1996 c 


NGC 3726 


H I 


Verheijen 1997 b 


NGC 3893 


Ha, H I 


Kranz 2002, Verheijen 1997 b 


NCC 3949 


H I 


Vorlnpiifin 1 QQ7 b 

V LI 11^1 I l_.ll _L C/ iJ I 


NGC 3953 


HI 


Verheijen 1997 b 


NGC 3992 


HI 


Verheijen 1997 b 


NGC 4051 


HI 


Verheijen 1997 b 


NGC 4062 


Ha 


Rubin ct al. 1985 


NGC 4138 


N II, H I 


Jore, Broeils, & Haynes 1996 c ,Verheijen 1997 b 


NGC 4651 


Ha 


Rubin, Waterman, & Kenney 1999 


NGC 4698 


Ha 


Rubin, Waterman, & Kenney 1999 


NGC 5371 


HI 


Bcgeman 1987 


NGC 5806 


Ha 


Courteau 1997 


NGC 6300 


FP Ha, H I 


Buta et al. 2001, Ryder et al. 1996 


NGC 7083 


Ha 


Mathewson, Ford, & Buchhorn 1992 


NGC 7217 


Ha 


Rubin et al. 1985 


NGC 7606 


Ha 


Mathewson, Ford, & Buchhorn 1992 



a The notation "FP Ha" is used for rotation curves derived from Fabry-Perot measurements of Ha. 
b Errors are taken as the difference in velocity between the approaching and receding sides. 
c The rotation curve has been extracted electronically from a plot in the referenced paper. 
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Table 2 





Observed 




Baryonic Matter 


Dark Matter 




Vtot ,max 


tT 

*b,max 


R{Vb,max) 






Rx 






/ \ i r\ i ^ ™ 
*+U. laex 






A+o.lderr 


A+o.lderr 






A n 1 fiprp 

' — 1 — U.J. cicx 




A-O.ldea; 


A-O.ldea; 


A-O.ldea; 






/\ n QWc-v 
*— » — u.oaea; 




A_o.3dex 


A-0.3dex 


A_o.3dex 


v_ < ( x i < ( ^A. y 




fkm / si 


(kpc) 


(1O 1O M ) 


(kpc) 


(kpc) 


NGC 157 


205 


205 


6.2 


7.1 










+25 




+ 1.8 










-22 




-1.4 










-60 




-3.5 






NGC 289 


182 


189 


3.6 


4.1 


9.1 


11.3 






+23 




+ 1.0 


+ 1.7 


+4.2 






-20 




-0.8 


-1.8 


-2.7 






-55 




-2.0 


-4.0 


-4.7 


NGC 488 


350 


296 


10.1 


31.1 










+36 




+7.4 










-32 




-5.9 










-86 




-14.3 






NGC 908 


200 


201 


7.3 


7.2 










+24 




+ 1.8 










-22 




-1.5 










-59 




-3.5 






NGC 1087 


136 


141 


5.9 


2.8 










+ 17 




+0.7 










-15 




-0.6 










-41 




-1.4 






NGC 1090 


170 


165 


7.3 


5.4 


16.4 


23.2 






+20 




+ 1.4 


+4.4 


+7.8 C 






-18 




-1.1 


-5.3 


-4.9 






-49 




-2.7 


-12.5 


-18.2 


NGC 1241 


300 


250 


7.9 


19.2 


12.4 








+30 




+4.8 


+4.4 d 








-27 




-3.8 


-5.6 








-73 




-9.2 


-12.4 d 




NGC 1385 


140 


133 


4.0 


2.6 










+ 16 




+0.6 










-14 




-0.5 










-39 




-1.2 






NGC 1559 


150 


136 


4.5 


2.3 










+ 17 




+0.6 










-15 




-0.5 










-40 




-1.1 
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Table 2 — Continued 





Observed 




Baryonic Matter 


Dark Matter 




'Tot TTLdX 


*c> : 7nax 


tt(Vb,max) 


M b 


-Rio 


i?X 






^+0.1dex 




A 

L±+0.1dex 


A 

L±+0.1dex 


A 

t-±+0.1dex 






A— o. ldex 




A-O.ldex 


A-O.ldex 


A-O.ldex 






^—0.3dex 




A_o.3dea; 


A_o.3de£c 


A-0.3dex 


Galaxy 


(km/s) 


(km/s) 


(kpc) 


(lO iU M ) 


(kpc) 


(kpc) 


NGC 1832 


200 


209 


3.8 


5.4 










+25 




+ 1 4 










-23 




-1.1 










-61 




-2.7 






NGC 2090 a 


160 


153 


1.6 


1.6 










+ 19 




-1-0 4 










-16 




-0.3 










-45 




-0.8 






NGC 2139 


140 


121 


4.5 


1.9 










+ 15 




+0.5 










-13 




-0.4 










-35 




-0.9 






NGC 2280 a 


210 


177 


6.9 


6.8 


8.7 








+5 




+ 1 8 


+ 1 8 








-19 




-1.4 


-3.5 








-32 




-3.4 


-5.4 




NGC 2841 


325 


284 


4.7 


13.0 


9.0 


12.6 






+49 




+3.3 


+3 4 


+5.8 






-20 




-2.6 


-0.6 


-2.9 






-35 




-6.4 


-6.2 


-3.7 


NGC 3198 


152 


120 


5.6 


2.3 


6.8 


10.5 






+ 14 






-to. 1 


4-9 fi 






-13 




-0.5 


-2.1 


-2.5 






-35 




-1.1 


-4.6 


-7.3 


NGC 3223 


320 


314 


10.4 


30.8 










+38 




+7.6 










-34 




-6.1 










-92 




-14.7 






NGC 3319 b 


132 


50 


9 


0.7 


2.6 


3.5 






+7 




+0.1 


+0.7 


+ 1.1 






-5 




-0.1 


-0.6 


-0.8 






-14 




-0.3 


-2.6 d 


-3.5 d 


NGC 3521 


221 


263 


3.1 


8.0 


10.3 


13.3 






+32 




+ 1.9 


+2.3 


+4.2 






-29 




-1.5 


-1.8 


-2.4 






-34 




-3.6 


-5.8 d 


-6.3 



34 



Table 2 — Continued 



Observed Baryonic Matter Dark Matter 





T 7 

vtot,max 


T T 

vb,max 


tt\Vb,max) 


71 J 




r> 

Kx 






A 




A 

Li+O.ldex 


A 

£±+0.1dex 


A 

L±+0.1dex 






A 

^-O.ldex 




A 

A-O.ldex 


A 

A-O.ldex 


A 

A—O.ldex 






A 

^-O.Sdex 




A 

A-0.3dex 


A 

A-0.3dex 


A 

A—0.3dex 


Galaxy 


(km/s) 


(km/s) 


(kpc) 


/1 nlO 71//" \ 

(1O 1U M ) 


(kpc) 


(kpc) 


NGC 3726 


169 


144 


10.8 


6.0 


20.4 


25.3 






+22 




+ 1.5 


+3.2 


+3.5 






-12 




-1.2 


-10.6 


-2.2 






-42 




-2.9 


-20.4 d 


-16.3 


NGC 3893 


210 


186 


5.3 


6.8 


13.3 


19.4 






+23 




+ 1.6 


+4.5 


+2.6 C 






-20 




-1.3 


-5.5 


-4.7 






-54 




-3.2 


-10.3 


-13.4 


NGC 3949 


169 


158 


3.9 


2.7 


7.4 








+ 18 




+0.6 


+0 7 d 








-17 




-0.5 


-1.6 








-46 




-1.2 


-7.4 d 




NGC 3953 


225 


227 


7.7 


12.0 


17.8 








+28 




+3.1 


+0.3 C 








-25 




-2.4 


-4.2 








-66 




-5.9 


-14.0 




NGC 3992 


272 


188 


13.1 


11.2 










+21 




+2.8 










-21 




-2.3 










-56 




-5.5 






NGC 4051 


170 


167 


0.3 


2.4 


5.8 








+20 




+0.6 


+2 1 








-18 




-0.5 


-1.3 








-49 




-1.1 


-5.8 d 




NGC 4062 


162 


110 


3.4 


1.1 


d 


3.9 






+ 14 




+0.3 


+ 1 


+2.4 






-11 




-0.2 


o d 


-2.0 






-32 




-0.6 


o d 


-3.9 d 


NGC 4138 


195 


272 


1.0 


4.2 


7.5 


17.8 






+33 




+ 1.0 


+8.9 


+3.1 






-30 




-0.8 


-1.4 


-5.8 






-79 




-2.0 


-7.5 d 


-12.6 


NGC 4651 


210 


173 


2.6 


3.7 










+32 




+0.9 










-10 




-0.8 










-45 




-1.8 
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Table 2 — Continued 





Observed 




Baryonic Matter 


Dark Matter 






^6, max 


Ri Vh 777 OT ) 




Rio 












ZA-un 1 dpT 

\J ■ J- L- ul> 


1\a-0 1 rtpT 


A 4-0 1 dpT 






A 

A-O.ldex 




A_ n i r/<=-r 


A_ n i dp-r 


A_ n i dp-r 






A_o.3de£c 




A— o 3dex 


A_o.3rfea: 


A_o.3rfea: 


Galaxy 


(km/s) 


(km/s) 


(kpc) 


(10 10 M Q ) 


(kpc) 


(kpc) 


NGC 4698 


220 


223 


3.8 


6.5 










+ 14 




+ 1.6 










-35 




-1.3 










-68 




-3.2 






NGC 5371 


242 


289 


13.9 


33.5 










+21 




+8.5 










-43 




-6.7 










-87 




-16.4 






NGC 5806 


200 


190 


1.1 


5.3 










+27 




+ 1.3 










-17 




-1.1 










-27 




-2.6 






NGC 6300 


208 


220 


4.3 


6.7 


13.7 


20.4 






+ 12 




+ 1.7 


i C O 

+5.8 


+2.6 






-35 




-1.4 


-1.3 


-5.5 






-67 




-3.3 


-8.8 


-13.5 


NGC 7083 


210 


245 


6.5 


13.1 










+30 




+3.3 










-27 




-2.7 










-72 




-6.4 






NGC 7217 


284 


283 


2.0 


7.9 










+35 




+2.0 










-31 




-1.6 










-84 




-3.9 






NGC 7606 


280 


228 


9.8 


12.9 


2.2 








+28 




+3.2 


+0.8 C 








-25 




-2.6 


-0.7 








-67 




-6.2 


-2.2 d 





a Quantities are calculated from B — V instead of B — R 
b Quantities are calculated from (M/L)*.h instead of {M/L)*_k 
c Upper limit 
d Lower limit 
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Table 3 





orig 


;onal color-M/L 
V200 x 2 




+0.1 dex 






-0.1 dcx 






-0.3 dex 






c 




^200 


y 2 




^200 


Y 2 
A 




V2Q0 


y 2 


Galaxy 




(km/s) 






(km/s) 






(km/s) 






(km/s) 




NGC 157 


0.5 


111 


9.5 


0.5 


61 a 


19.9 


4.0 


96 a 


6.7 


16.0 


96 


5.8 




0.5 


61 a 


13.8 


0.5 


61 a 


28.2 


0.5 


61 a 


8.2 


2.0 


131 


7.3 


NGC 289 


8.0 


143 


0.7 


6.5 


148 


0.8 


10.5 


138 


0.8 


11.5 


138 


0.8 




3.0 


168 


1.2 


1.0 


238 


1.7 


5.0 


133 


0.8 


7.5 


143 


0.7 


NGC 488 


16.5 


199 a 


9.1 


8.0 


204 


6.5 


20.0 


214 a 


12.6 


20.0 


279 


22.7 




3.0 


204 


5.3 


0.5 


199 a 


6.8 


8.0 


199 a 


5.2 


18.5 


199 a 


6.1 


NGC 908 




















20.0 
14.0 


116 

96 a 


8.7 
6.5 


NGC 1087 




















16.5 
4.0 


68 a 
68 a 


1.3 
1.5 


NGC 1090 


1.0 


400 


6.4 


0.5 


400 


9.7 


5.0 


160 


5.3 


12.0 


125 


4.7 




0.5 


105 


11.4 


0.5 


80 a 


15.6 


0.5 


215 


8.7 


1.0 


360 


6.6 


NGC 1241 


1.5 


745 


4.7 


0.5 


150 a 


9.3 


3.0 


745 


3.4 


20.0 


180 


3.4 




0.5 


150 a 


7.3 


0.5 


150 a 


22.2 


0.5 


354 


4.0 


1.5 


585 


2.3 


NGC 1385 




















4.0 
0.5 


328 
143 


0.1 
0.7 


NGC 1559 




















3.5 
0.5 


350 a 
350 a 


5.1 
7.0 


NGC 1832 


0.5 


89 a 


8.8 








0.5 


394 


4.2 


20.0 


94 


4.0 




0.5 


89 a 


19.1 








0.5 


394 


6.9 


0.5 


389 a 


3.7 


NGC 2090 














20.0 
2.5 


77 a 
77 a 


0.8 
0.4 


20.0 
13.5 


122 
77 a 


1.2 
0.3 


NGC 2139 


2.5 


301 


1.5 


0.5 


71 a 


2.2 


4.0 


336 


1.1 


7.5 


261 


0.8 




0.5 




2.7 


0.5 


71 a 


5.6 


0.5 


156 


2.5 


2.0 


336 


2.3 


NGC 2280 


2.0 


470 


1.2 


1.5 


475 


0.5 


3.0 


495 


0.2 


5.0 


375 


0.2 




0.5 


105 a 


3.9 


0.5 


105 a 


2.4 


0.5 


265 


2.7 


1.0 


525 


2.6 


NGC 2841 


8.5 


242 


13.4 


3.0 


402 


27.6 


12.5 


217 


11.0 


15.0 


207 


7.0 




0.5 


582 


23.3 


0.5 


452 a 


53.9 


1.0 


647 


18.6 


3.5 


337 


11.6 


NGC 3198 


6.0 


140 


1.7 


4.0 


160 


2.8 


8.0 


130 


1.1 


11.5 


120 


0.6 




0.5 


345 


4.0 


0.5 


285 


6.5 


1.5 


235 


2.8 


5.0 


145 


1.5 


NGC 3223 


0.5 


733 


5.3 


0.5 


158 a 


14.8 


9.5 


178 


5.0 


20.0 


178 


7.1 




0.5 


158 a 


10.1 


0.5 


158 a 


39.9 


0.5 


233 


3.7 


11.5 


158 


2.6 


NGC 3319 


2.0 


269 


0.3 


1.5 


314 


0.3 


2.5 


233 


0.2 


3.5 


189 


0.2 




1.0 


304 


0.9 


1.5 


314 


1.3 


1.5 


259 


0.7 


1.5 


294 


0.4 


NGC 3521 


12.5 


121 


3.1 


7.0 


136 


4.3 


18.5 


116 


2.6 


20.0 


121 


2.2 




0.5 


281 


4.1 


0.5 


171 


6.2 


4.5 


146 


3.5 


13.5 


116 


2.1 


NGC 3726 


1.0 


334 


1.0 


0.5 


379 


3.6 


1.0 


394 


0.3 


6.0 


154 


0.6 




0.5 


159 


5.9 


0.5 


94 


10.2 


0.5 


239 


2.9 


0.5 


414 


0.4 
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Table 3 — Continued 



Galaxy 


orig 

c 


;onal color-M/L 
V 2 qo x 2 
(km/s) 


c 


+0.1 dex 

V200 
(km/s) 


x 2 


c 


-0.1 dex 

V200 
(km/s) 


x 2 


c 


-0.3 dex 

V200 
(km/s) 


X 2 


NGC 3893 


4.0 


159 


2.3 


0.5 


334 


4.4 


18.0 


94 


2.2 


20.0 


119 


4.0 




0.5 


74 


3.7 


0.5 


74a 


8.9 


0.5 


204 


2.3 


8.5 


119 


2.1 


NGC 3949 


2.0 


390 


2.2 


0.5 


365 


5.5 


3.0 


405 


0.7 


7.0 


245 


0.03 




0.5 


85 a 


4.5 


0.5 


85 a 


10.8 


0.5 


280 


2.9 


2.0 


385 


0.8 


NGC 3953 


0.5 


563 a 


2.0 


0.5 


113 a 


6.5 


1.5 


558 


0.3 


9.5 


183 


0.02 




0.5 


113 a 


5.3 


0.5 


113 a 


22.9 


1.5 


183 


2.2 


1.0 


493 


0.2 


NGC 3992 


20.0 


149 


0.3 


14.0 


154 


0.3 


19.5 


159 


0.4 


20.0 


164 


1.2 




8.0 


164 


0.3 


2.0 


234 


0.3 


11.0 


164 


0.3 


20.0 


153 


0.4 


NGC 4051 


2.0 


382 


1.9 


1.0 


322 


2.9 


3.0 


352 


1.3 


10.0 


152 


1.0 




0.5 


137 


3.8 


0.5 




4.9 


0.5 


322 


2.5 


2.0 


352 


1.7 


NGC 4062 


18.5 


97 


1.7 


3.5 


382 


1.5 


20.0 


102 


1.7 


20.0 


117 


2.0 




1.5 


372 


0.6 


0.5 


352 a 


1.1 


2.5 


357 


0.5 


4.0 


377 


0.7 


NGC 4138 


1.0 


375 a 


15.7 


0.5 


295 


30.1 


8.5 


115 


9.9 


20.0 


105 


7.3 




0.5 


75 a 


23.6 


0.5 


75 a 


43.0 


0.5 


75 a 


15.4 


1.0 


375 


10.8 


NGC 4651 


3.0 


517 


2.3 


1.0 


557 a 


4.9 


15.0 


142 


2.0 


20.0 


152 


3.0 




0.5 


112 a 


6.1 


0.5 


112 a 


11.1 


1.0 


307 


5.8 


3.5 


382 


5.7 






















O.O 


^9^ 
DZO 


1 Q 9 

iy .z 






















0.5 


105 a 


29.2 


NGC 5371 














0.5 


217 


6.2 


9.0 


137 


0.9 
















0.5 


107 a 


15.1 


0.5 


197 


0.8 


NGC 5806 














14.0 


96 a 


1.6 


20.0 


96 a 


3.7 
















1.0 


101 


0.7 


6.0 


96 a 


1.2 


NGC 6300 














5.0 


174 


9.5 


20.0 


109 


9.3 
















0.5 


209 a 


7.0 


8.0 


119 a 


4.6 


NGC 7083 




















13.0 


126 


1.1 






















0.5 


361 


0.9 


NGC 7217 




















20.0 


176 


15.6 






















19.0 


106 a 


4.0 


NGC 7606 


17.5 


142 a 


1.7 


9.5 


142 a 


2.5 


20.0 


157 


1.4 


20.0 


192 


2.0 




0.5 


447 


3.0 


0.5 


142 a 


3.6 


2.5 


357 


2.7 


12.0 


187 


2.3 



a Limit of parameters searched. 



Note. — The first and second rows for each galaxy list information for the uncontracted and contracted 
halo models. 
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Table 4 







Correlation Coefficients 




Galaxy Property 


All Data 


Without Outlicrs a 




Rx/hiR Rx 


Rx/hm 


Rx 


vtot^max 


0.41 




0.90 




Vb,max 


0.85 




0.96 




R{Vb,max) 




0.22 




0.58 b 


T-type 


0.89 




0.90 




logi M b 


0.65 




0.90 








0.21 




0.80 c 


Rw 




0.94 




0.98 


B 


0.36 




0.77 




K d 


0.41 




0.80 






0.51 




0.68 




Mo, if 


0.64 




0.87 





a Outliers are NGC 2841 and NGC 4138. 

b Outliers for this relation arc NGC 3319 and NGC 4062. 

C NGC3319 is not included because it is an outlier in this relation. 

d NGC 3319 is not included because it does not have if -band surface brightness profiles with high enough 
signal-to-noise to derive an integrated magnitude. 
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